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This  report  presents  and  analyzes  simulated  results  for  electric  field  propa- 
gation through  a thick,  extended  region  of  high-altitude  striationi  electron! 
density.  The  simulated  results  correspond  to  the  problem  of  communication  fronj 
a )SCS  satellite  to  a ground  station  at  a frequency  of  7.5  GHz.  Statistical 
re>ults  are  presented  for  S4  scintillation  index  and  intensity  correlation 
length  for  a power- law  and  an  exponential  phase-screen  PSD  for  values  of  phase 
standard,  deviation  arranging  from  1 to  100  radians  and  phase  correlation 
length  -pranging  from  1 to  10  km.  \ 
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20.  Abstract  (Continued) 


Examples  of  specific  realizations  of  the  electric  field  amplitude  and  phase 
at  the  receiver  are  also  shown.  These  realizations  of  signal  structure  are 
being  used  in  simulations  of  DSCS  receiver  operation,  reported  elsewhere. 
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SECTION  1 
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1.0  INTRODUCTION 

The  problem  of  communication  from  a satellite  to  a ground  station 
in  the  presence  of  large,  spatially  extended  regions  of  high-altitude  nuclear 
burst  produced  striations  of  electron  density  can  be  modeled  as  the  problem 
of  propagation  through  a thick  medium  composed  of  random  index-of-refraction 
fluctuations.  Since  no  general  analytical  solution  is  available  for  this 
problem  it  must  be  landled  numerically.  In  this  report,  a multiple  phase- 
screen  propagation  simulation  is  developed  as  a method  to  propagate  electro- 
magnetic waves  through  statistically  chosen  realizations  of  the  random  medium. 

In  general,  the  electron  density  striations  are  represented  by  a 
luunber  of  phase  screens  which  are  chosen  to  possess  the  same  stati.stical 
description  as  tl>e  original  striations.  A i)lane  wave  is  then  propagated 
from  one  phase  screen  to  the  next  until  a solution  for  the  electric  field 
in  the  receiver  plane  is  obtained.  To  obtain  the  statistics  of  the  received 
field,  the  problem  is  solved  a number  of  times  with  different  realizations 
of  the  phase  screens  (all  with  the  same  statistics)  and  the  results  are 
averaged. 


The  possibility  of  generating  and  retaining  the  realizations  of 
the  received  electric  field  is  an  advantage  of  this  type  of  simulation 
solution  particularly  when  the  i-esults  are  to  be  used  in  conjunction  with 
a receiver  simulation  of  the  ground  station  electronics.  For  all  the  work 
reported  here  several  realizations  of  the  received  field  were  retained  and 
are  available  for  future  analysis. 
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This  report  considers  the  effect  of  nuclear  burst  produced  stri- 
ation  structure  on  propagation  of  7.5  GHz  (wavelength  = 0.04  in)  communica- 
tion signals  from  a satellite  to  a ground  station.  Die  two  different 
geometries  relevant  to  the  Defense  Satellite  Communications  System  (DSCS) 
considered  are  shown  in  Figure  1 (a)  for  the  case  of  striations  ranging 
in  altitude  from  300  km  (F  region  altitude)  to  15000  km  and  ca  e (b) 
for  striations  ranging  from  6000  km  to  15000  km  altitude,  fn  both  cases  the 
striations  are  assumed  to  be  elongated  in  the  direction  perpendicular  to 
the  page  so  that  there  is  no  variation  in  the  y-direction.  All  results  are 
given  in  terms  of  the  one  dimensional  phase-screen  power  spectral  density  or 
autocorrelation  function.  Tlie  relationship  of  these  quantities  to  the  elec- 
tron density  power  spectral  density  is  given  in  Appendix  B.  The  choice  of 
phase-screen  power  spectral  density  is  arbitrary  but  for  this  work  results 
were  obtained  for  two  spectra: 

(...) 

- (i.r  i,p 


f 9 \ ">  *‘0  “ ! ^ I ^0 


As  shown  in  Appendix  B ^ corresponds  to  a K " power  law  spectral 

density  for  one-dimensional  electron  density  fluctuations  (thought  to  be 

(2) 

representative  of  the  ambient  ionosphere)  and  corresponds  to  a one- 

dimensional exponential  (or  Chesnut)  power  spectral  density  for  integrated 
electron  density  fluctuations. 

The  corresponding  phase  screen  autocori'elation  functions  may  be 
anal>tically  obtained  as: 


(1-3) 


B^^C) 


•f-H) 
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^ Satellite 
Transmitter 


Altitude  = 15000  km 


300  km 


Receiver 


Altitude  = 15000  km 


6000  km 


Receiver 


(a)  Striations  located  between 
300  km  and  15000  kn  altitude 


(b)  Striations  located  between 
6000  km  and  15000  km  altitude 


Figure  1.  Geometry  of  DSCS  simulation. 
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where  Kj^(x)  is  the  modified  Bessel  function  (Abramowitz  and  Stegun,  p.  374). 
(The  Fourier  transform  relationship  used  here  between  autocorrelation  func- 
tion and  power  spectral  density  is  given  by  Equations  B-11  and  B-14.)  For 
this  report  an  important  parameter  of  interest  is  the  autocorrelation  length 
which  is  defined  as  the  e ^ point  of  the  autocorrelation  function.  For 
the  power  law  spectrum,  4*^^^  the  coi'relation  length  is  = 1.55  L^, 

and  for  the  exponential  spectrum,  the  correlation  length  is  given 

(2) 

by  - i*31  Lq.  For  this  report,  the  correlation  length  ranged  from 

1 to  10  kilometers  for  both  choices  of  the  phase-screen  power  spectral 
density. 


The  final  physical  parameter  of  interest  is  the  total  phase  stan 
dard  deviation  of  the  phase-screens  comprising  the  multiple  phase-screen. 
Since  each  phase-screen  is  designed  as  an  assembly  of  Gaussian  .'andom  vari 

7 

ables  with  zero  moan  and  variance  a~  and  is  independent  of  the  other 
individual  phase-screens,  the  total  variance  is  the  sum  of  the  individual 
variances  and  thus  the  phase  standard  deviation  is  given  by 


(1-5) 


where  N is  the  total  number  of  phase-screens.  For  this  work  N was 
chosen  to  be  10  and  ranged  from  l.O  to  100  radians. 

4) 


The  main  body  of  tins  report  is  divided  into  three  major  sections. 
Sections  2 and  3 show  the  application  of  the  MPS  code  to  problems  whose 
analytic  solutions  are  available  and  allow  a comparison  cf  the  code  results 
to  analytical  results.  Section  2 deals  with  the  calculation  of  the  dif- 
fraction pattern  of  a single  Gaussian  lens  in  the  center  of  a phase  screen. 
The  analytic  solution  is  obtained  as  a Fresnel-Kirchhoff  integral  which  may 
be  written  as  an  infinite  series  and  summed.  Section  3 shows  the  results 
for  a single  phase-screen  characterized  by  a Gaussian  power  spectral  density 
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and  small  phase  variance  a^.  Results  for  the  observed  intensity  power 
spectral  density,  and  scintillation  index  are  compared  with  analytic  cal- 
culations in  the  weak-scatter  approximation.  Section  4 describes  the 
results  of  the  bSCS  simulation  at  a frequency  of  7,5  GHz,  Appendix  A 
contains  a description  of  the  multiple  phase-screen  code,  its  analytical 
foundation  and  the  technique  used  to  generate  the  random  phase-screens. 
Appendix  B contains  a derivation  of  the  relationship  between  the  phase- 
screen  standard  deviation  o.  and  the  electron  density  fluctuations  and 
a brief  discussion  of  the  relationship  between  one,  two  and  three-dimensional 
power  spectra.  Appendix  C shows  one  realization  of  all  ten  phase-screens 
and  traces  the  signal  intensity  as  the  wave  propagates  through  the  medium 
for  a particular  geometry  relevant  to  DSCS. 
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SlCTION  2 

2.0  DIFFRACTION  BY  A GAUSSIAN  LENS 

2.1  Gaussian  Lens  - Code  Calculation 

As  a test  of  the  multiple  phase-screen  propagation  code,  a single 
Qiussian  phase  lens  given  by 

<>(x)  = exp(-x'/r')  (2-1) 

was  used  at  the  first  phase-screen  location  (z  = 0.0)w*th  4'.^  chosen  as 
10  radians  and  r^  taken  as  X. 

The  intensity  I = |I;|  was  obsorveu  as  a function  of  distance 
from  tlie  lens,  z/X  as  indicated  in  I'igurc  J(a-b),  !n  this  figure,  all 
tl>e  pluise -screens  and  ol)sorvation  screens  are  aligned,  so  that  features 
of  the  diffraction  pattern  can  be  easily  followed  for  changing  z/X  values. 

The  focal  length  for  a tJaussian  lens  is  given  by  Salpeter  (1907)  as 
kr“/2ip^  where  k is  the  wavenumber  2ii/X.  Taking  = X and  = 10  radians, 
the  focal  length  I'/X  is  0.31  which  corresponds  reasonably  well  with  the 
value  of  z/X  = 0.5  where  the  intensity  in  the  diffraction  pattern  at 
X = 0 builds  up  to  a maximum.  At  values  of  z/X  greater  than  the  focal 
length,  the  diffraction  pattern  exhibits  increasingly  more  complex  patterns 
associated  with  rays  coming  from  the  edges  of  the  lens,  rather  than  from 
the  center. 
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2.2 


Gaussian  Lens  - Analytic  Results 


An  analytical  relationship  between  the  electric  field  at  z = z^ 
and  z = z^  is  given  by  the  Frcsnel-Kirchhoff  integral  (Ratcliffe,  1956) 


- 1/9  CO 

i27r(z2-zp^  “ ^ -ikCz^-zp 


...... /. 


X exp  I i ~ k(x-?)“/(z2-Zj)  | n(C,zp 


(2-2) 


where 


l;(x,Zj)  is  the  electric  field  as  a function  of  x in  the  z = z^ 
plane.  I'or  an  initial  electric  field  in  the  z = 0 plane  given  by 

9 0 

(i(C  ) = exp  i {<{',,  expC-C^/r!!) },  l:t|uation  2-2  may  be  written  as  an  infini‘  ! 
u v)  9 2 

series  by  expanding  exp  i {4>q  cxp(-C“/r^))  in  a Taylor  series  as 


I 


n=0 


a.' 


(2-3) 


The  resulting  integral  over  may  then  be  analytically  performed 
using  liquation  3.323(2)  on  page  307  of  (Iradshteyn  and  Ryzhik  (1965).  The 
result  may  be  expressed  as 

i2'nz'  ,2  , (i‘l>n)'^  o i /-> 

h(x',z')  = /n  e 2-  ^ 


exp 


nz' 


-TI  x'^/z 


n=0 


n ! 


,2 


(n/r^‘'-iTr/z') 


(2-4) 


where  the  primed  quantities  x'  = v/A,  z'  = z/A,  and  r^  = r^/A,  This 
series  is  easily  summed  (numerically)  and  the  results  for  several  values 
of  z/A  are  shown  in  Figure  2 to  be  identical  to  the  code  calculation  for 
z/A  ^ 2.0.  11 


Phase  Pattern  of 
Gaussian  Lens 


Numerical  Results 


Theoretical  Results 


Distance  (x/A)  Distance  (x/A) 

Figure  2(a).  Diffraction  pattern  of  a Gaussian  lens  - numerical  results 
from  simulation  on  left-hand  side  versus  theoretical  results 
on  right-hand  side. 
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Theoretical  Results 


-10  -5  0 5 10 


Distance  (x/A) 


Distance  (x/A) 


Figure  2(b).  Diffraction  pattern  of  a Gaussian  lens  - numerical  results 
from  simulation  on  left-hand  side  versus  theoretical  results 
on  right -hand  side. 
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For  z/X  = 10  the  multiple  phase-screen  calculations  deviate 
from  the  theoretical  results  because  energy  which  has  left  one  side  of  the 
grid  is  coming  back  into  the  other  side  (this  is  the  well  known  wrap-around 
or  aliasing  effect  in  fast  Fourier  transforms).  Hendrick  (1977)  has  deter- 
mined criteria  required  to  avoid  this  and  other  numerical  problems  and  these 
criteria  were  used  to  assure  proper  application  of  the  code  results. 


SECTION  3 


3.0 


DIFFRACTION  BY  A RANDOM  THIN  PHASE-SCREEN 


This  section  contains  a second  application  of  the  multiple  phase- 
sci’een  propagation  code  to  a problem  for  which  analytical  results  are  avail- 
able - scattering  through  a thin  phase-screen  charactoriiied  by  a Gaussian 
power  spectral  density.  For  this  example,  only  one  phase-screen  was  used 
to  characterize  the  random  medium.  Tiie  screen  was  located  at  z = 0 and 
was  generated  with  a phase  power  spectral  density  (PSD)  of 

(3-1) 


where  the  phase  standard  deviation  was  chosen  as  = 0.1  radians  and  was 
taken  as  U,  The  corresponding  phase  autocorrelation  function  is  then 


B^(C)  = exp(-r/x^) 


(3-2) 


To  obtain  statistical  results  for  the  received  signal,  the  simu- 
lation was  exercised  ten  times  with  ten  different  random  phase-screens 
(each  screen  chosen  used  a different  set  of  random  numbers;  see  Appendix  A). 
Figure  3 shows  a comparison  of  the  intended  phase-screen  power  spectral 
density  where  K'  = KA  (dotted  line)  with  the  actual  mean  power 

spectral  density  obtained  by  averaging  the  power  spectral  densities  of  each 
of  the  ten  realizations.  Although  no  statistical  study  was  performed  to 
show  the  effect  of  averaging  over  different  numbers  of  realizations  this 
result  indicates  that  ten  realizations  are  sufficient  to  obtain  reasonably 
accurate  statistics. 


Intended  PSD 
(dashed  curve) 


i 10‘ 


--12F 


h!crma1 


Figure  3.  Comparison  of  intended  phai 
that  generated  by  the  simu 


2 

Figure  4 shows  the  average  intensity  (|e|  ) power  spectral  d^,Tisity 
at  the  locations  z/X  = 0.628,  6.28,  and  37.7,  each  compared  to  the  theoreti- 
cal results  from  Salpeter  (1967) 

4>j(K)  = 4<I>^(K)  sin^K^;;/2k)).  (3-3) 

(Strictly  speaking  is  the  PSD  of  the  normalized  intensity  deviation 
(|e1  - e| '^))/^ I e1  ^^) . Excellent  agreement  between  the  theoretical  and  simu- 
lation results  is  shown  with  the  simulation  results  for  z/X  = 6.28  even 
matching  the  location  of  the  nulls  (but  not  the  infinite  depth)  of  the 
analytical  result.  In  all  three  curves  the  simulation  deviates  from  the 
analytical  results  for  large  wavenumbers,  which  is  likely  caused  by  aliasing 
which  always  occurs  in  the  application  of  fast  Fourier  transform  techniques, 

4 second  possible  explanation  for  this  deviation  at  large  wavenumbers  is  that 
the  small  scale  sizes  may  produce  angular  scatter  which  causes  this  portion 
of  the  wave  to  alias  to  the  other  side  of  the  phase  screen.  However,  this 
aliasing  behavior  should  become  more  pronounced  as  the  pi'opagation  distance 
z/X  increases.  Since  this  is  not  the  case  for  Figure  4,  this  explanation 
must  be  discarded.  For  z/X  = 37.7  the  simulation  accurately  models  the 
peaks  of  the  analytical  curve.  However,  except  for  the  iirst  null  at 
KX  ~ 1 (K  z/2k  = ir),  the  details  of  the  nulls  are  lacking.  This  smearing 
of  the  nulls  is  caused  by  the  averaging  process. 

The  circled  dots  of  Figure  5 show  the  value  of  scintillation 
index  obtained  from  the  simulation  for  values  of  of  0.1,  1,0  and 
10,0  radians.  The  solid  curves  are  analytical  estimates  taken  from  the 
weak-scatter  approximation  of  Salpeter 
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Figure  5.  Circled  dots  show  S4  scintillation  index  for  Gaussian 
spectrum.  The  solid  lines  are  analytic  approximations 
valid  for  the  weak-scatter  case  {Xq  = X). 


which  may  be  exactly  integrated  for  the  Gaussian  phase  power  spectral 
density  given  by  Equation  (3-1)  to  obtain 


It  is  seen  from  Figure  5 that  excellent  agreement  is  obtained  between  this 
analytical  approximation  and  simulation  results  when  the  weak-scatter 
assumption  inlierent  in  Equation  3-5  is  valid. 


20 


SECTION  4 


4.0  DSCS  SIMULATION  RESULTS 

This  section  describes  the  results  of  the  simulation  as  applied  to 
the  satellite  communication  problem  discussed  in  the  introduction.  As  stated 
there,  the  simulation  was  applied  to  the  two  geometries  shown  in  Figure  1, 
for  two  different  phase-screen  PSD's  as  given  by  Equations  1-1  ariU  1-2.  The 
phase  correlation  lengths  ranged  from  1 to  10  km,  and  the  phase  standard 
deviation  (as  given  by  Equation  1-5)  ranged  from  1.0  to  100  radians.  Table 
4-1 (a)  and  (b)  shows  the  cases  which  were  simulated  for  the  power-law  PSD 
and  for  the  exponential  PSD  respectively.  In  addition  to  the  phase  correla- 
tion length,  It^,  and  phase  standai'd  deviation,  o^,  for  each  case  the  simu- 
lated pljaso-screen  length  L is  also  included.  In  all  cases,  for  this 
section  of  the  report,  the  statistical  results  were  obtained  by  averaging 
results  for  ten  different  realizations  of  the  problem.  For  each  realiza- 
tion, the  striated  region  was  I'epresented  by  ten  phase-screens,  each  com- 
posed of  2048  points. 

In  general,  the  results  were  not  a strong  function  of  the  two 
geometries  shown  in  Figure  1.  For  the  weak  scattering  cases  listed  in 
Table  4,  the  scintillation  index  for  the  thicker  striation  region  geometry 
(Figure  1(a))  was  up  to  25%  less  than  that  for  the  thinner  region  with  the 
same  total  phase  standard  deviation.  This  is  intuitively  expected  since 
for  the  thicker  geometry,  the  phase-screens  very  close  to  the  receiver 
plane  will  not  be  as  effective  at  producing  amplitude  scintillation  as 
phase-screens  farther  away.  In  the  rest  of  this  section  only  results  for 
the  geometry  of  Figure  1(a)  are  presented. 
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Table  4-1 (a).  Simulation  parameters  for  power-lav;  PSD. 


Phase  Correlation  Phase  Phase  Screen  Resulting 

Length  Standard  Deviation  Length  Scintillation  Index 

Z a.  I Geometry  1*  Geometry  2** 

c 2 


1 km 

1.0  rad 

30  km 

.37 

.42 

1 

2.5 

30 

.83 

.85 

1 

3.16 

30 

.94 

.93 

1 

10.0 

30 

1.05 

.98 

3 

3.16 

50 

.41 

.48 

3 

10.0 

50 

1.03 

1.06 

10 

10.0 

100 

.41 

.49 

10 

31.6 

100 

1.05 

1.05 

Table  4-1 (b).  Simulation 

parameters  for 

exponential  PSD. 

Phase  Correlation 

Phase 

Phase  Screen  Resulting 

Length 

Standard  Deviation  Length 

Scintillation 

Index 

^c 

2 

L 

Geometry  1* 

Geometry  2** 

1 km 

1.0  rad 

30  km 

.32 

.39 

1 

2.5 

30 

.79 

.87 

1 

3.16 

30 

.93 

.96 

1 

10.0 

30 

1.15 

1.04 

3 

3.16 

50 

.16 

.20 

3 

10.0 

50 

.67 

.89 

3 

31.6 

50 

1.60 

1.48 

10 

10.0 

100 

.05 

.06 

10 

31.6 

100 

.15 

.19 

10 

50.0 

100 

.25 

.34 

10 

100.0 

100 

.89 

1.15 

Geometry  1 refers  to  Figure  1(a) 
Geometry  2 refers  to  Figure  1(b). 
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In  Figure  6 the  scintillation  index  is  shown  as  a function  of  the 
total  integrated  phase  standard  deviation  for  both  spectra  for  values 
of  the  striation  autocorrelation  distance  of  1,  3 and  10  kilometers.  As 
expected,  the  scintillation  index  increases  with  increasing  0^  but  should 
ultimately  saturate  at  a value  of  unity,  characteristic  of  Rayleigh  electric 
field  strength  statistics  (Fante,  1977).  Values  greater  than  unity  prior  to 
saturation  are  not  prohibited  and  have  been  observed  in  ionospheric  measure- 
ments (l)v.  n.J.  Fremouw,  personal  conununication,  Stanford  Research  Institute,  1977) 
The  large  value  of  1.6  for  Sj  in  Figure  6 is  caused  by  strong  focussing  effects. 

Figure  7 shows  results  obtained  for  the  autocorrelation  distance 
of  the  signal  intensity  as  a function  of  the  total  phase  standard  deviation 
for  the  same  cases  as  shown  in  Figure  6.  These  results  wore  obtained  by 
finding  the  distance  corresponding  to  the  o~^  point  from  the  circular 
correlation  functions  obtained  from  the  average  of  the  ten  realizations. 

The  general  trend  of  decreasing  correlation  distance  with  increasing 
is  expected  intuitively  as  is  the  association  of  large  autocorrelation 
length  with  large  striation  autocorrelation  distance.  The  range  of  these 
curves  is  limited  bv  roquii’cments  on  the  phase-screen  sampling  distances. 

In  order  to  extend  the  calculation  of  the  intensity  correlation  distance 
farther  into  the  multiple-scatter  I'cgirae  two  techniques  appear  useful. 

The  first  involves  a numerical  solution  of  the  differential  equations  for 
the  mutual  coherence  function  (l(p)  I(p+C))  as  discussed  by  Yeh,  et  al., 

(1975).  A second  possible  method  involves  an  approximate  theory  developed 
by  Fante  (1975).  Both  techniques  require  further  evaluation. 
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Scintillation  Index 
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Phase  Standard  Deviation  o.  (radians) 


Figure  7.  Signal  intensity  correlation  distance  versus  as  a 
function  of  phase-screen  PSD  and  correlation  length  l 
for  «._  = !,  3,  and  10  km. 


4.1 


Simulation  Results  - 1 kilometer  Correlation  Length 


This  section  and  the  following  two  sections  discuss  the  results 
obtained  for  realizations  of  the  DSCS  propagation  simulation.  In  all  cases 
the  random  number  seed  was  the  same  for  eacl\  computer  run  so  that  the  set 
of  phase-screen  realizations  was  similar  in  "overall"  appearance  for  all  the 
runs.  This  also  has  the  result  that  some  features  of  the  received  intensity 
and  phase  are  recognizable  as  the  phase  standard  deviation  and  correlation 
length  change. 

Figure  8 shows  a comparison  of  the  actual  power  spectral  densities 
of  the  phase-screens  as  generated  by  the  MPS  simulation  for  a one  kilometer 
correlation  length.  Since  the  Fresnel  length  /X?  ranges  from  0.1  km 
(z  = 300  km)  to  0.8  km  (z  = 15000  km)  which  correspond  to  the  wavenumber 
regime  from  7.8  km"^  to  62.8  km'^  whei’o  the  spectra  are  somewhat  similar, 
the  overall  simulation  results  should  be  similar  for  those  two  spectra, 
since  amplitude  scintillations  are  caused  by  irregularities  whose  scale 
sizes  are  of  the  order  of  the  Fresnel  length.  But  since  the  power- law 
spectrum  is  somewhat  higher  at  high  wavenumbors  (corresponding  to  smaller 
scales),  the  received  electric  field  should  exhibit  more  small-scale 
structure  for  the  power- law  PSD. 

This  behavior  is  evident  in  Figures  9 and  10  which  show  one 

2 

realization  of  received  intensity  (I  = |n|  ) for  the  two  PSD's  for 
values  of  the  phase  standai-d  deviation  of  1.0,  3,16  and  10.0  radians. 

Figures  11  and  12  show  the  phase  (plotted  between  ±ir)  coi'responding  to  the 
intensities  shown  in  Figures  9 and  10. 

Figures  13  and  14  show  the  intensity  power  spectral  density  cor- 
responding to  the  cases  shown  in  Figures  9 and  10.  (The  results  shown  in 
Figures  13  and  14  result  from  averaging  the  PSD's  of  ten  realizations,  the 
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Figure  8.  Comparison  of  simulation  generated  phase-screen  PSD' 
correlation  length  = 1 km,  A = 0.04  m.  {a(K.  = /nj 
radians  in  figure).  ’ 
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Figure  9.  Realizations  of  signal  intensity  for  power-law 
PSD  with  = 1 km. 
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Figure  11.  Realizations  of  signal  phase  for  power-law  PSD 
with  = 1 km. 
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Figure  12.  Realizations  of  signal  phase  for  exponential  PSD  with 
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Normalized  Intensity  PSD  <{),(K)/X 


Figure  14.  Normalized  intensity  PSD  for  exponential  phase-screen 
PSD  with  " 1 km,  X = 0.04  m. 
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first  of  which  is  shown  in  Figures  9 and  10.)  For  the  weak  scattering  cases 
(o^  = 1.0)  the  intensity  PSD  closely  matches  the  PSD  of  the  phase-screen 
for  wavenumbers  greater  than  the  inverse  Fresnel  length.  This  behavior  is 
required  by  Equation  3-3  for  weak  scattering  since  the  intensity  PSD  is 
the  product  of  the  phase-screen  PSD  with  a sine-squared  factor  which  averages 
to  1/2  for  wavenumbers  much  greater  than  the  inverse  Fresnel  length.  For 
the  strong  scattering  cases  (o^  > 1.0),  the  intensity  PSD  gradually  becomes 
flatter  near  the  inverse  Fresnel  length  but  retains  a high-wavenumber  roll- 
off where  the  spectrum  again  has  the  same  slope  as  the  phase-screen  PSD. 

This  flat  low-wavenumber  spectrum  with  a high-wavenumber  rolloff  is  similar 
to  observations  of  intense  ionosphei'ic  scintillations  shown  by  Whitney  and 
Basu  (1977),  This  flattening  of  the  intensity  PSD,  of  course,  also  coi'res- 
ponds  to  a sharpening  of  the  intensity  autocorrelation  function  with  a 
resulting  decrease  in  the  signal  intensity  correlation  distance  with  in- 
creasing 0^.  This  behavior  is  shown  in  Figure  7. 
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4.2 


Simulation  Results  - 3 kilometer  Correlation  Length 


Figui'C  15  shows  a comparison  of  the  measured  (simulation)  phase- 

screen  PSD's  for  a phase  correlation  length  of  3 km.  For  this  case,  the 

power- law  PSD  is  much  larger  than  the  exponential  PSD  near  wavenumbers 

corresponding  to  the  Fresnel  length,  and  therefore  one  would  expect  a mcch 

greater  difference  between  the  received  signals  for  the  = 3 km  case 

than  for  the  £ = 1 km  case  discussed  in  the  previous  section, 
c 

The  difference  is  apparent  in  Figures  16  and  17  which  show  realiza- 
tions of  the  received  intensity  for  values  of  = 3.16  and  10.0  radians 
for  the  power-law  PSD  and  = 3.16,  10.0  and  31.6  I'adians  for  the  exponen- 
tial PSD.  Note  the  absence  of  small  scale  structure  for  the  exponential  PSD 

for  o.  = 3.16  and  10.0  radians. 

(I* 

The  phases  cori'osponding  to  the  intensities  of  Figures  16  and  17 
are  shown  in  Figures  18  and  19. 

For  the  interested  reader  Appendix  C shows  the  realization  of  the 
ten  phase-screens  and  the  electric  field  intensity  as  the  wave  progresses 
through  the  screens  to  the  observation  plane  for  the  case  = 31,6  and 
the  exponential  PSD. 

Figures  20  and  21  show  the  mean  intensity  PSD’s  corresponding  to 
the  intensity  realizations  shown  in  Figures  16  and  17.  Again,  as  the  phase 
standard  deviation  increases,  the  intensity  PSD  goes  from  a close 
resemblance  to  the  phase-screen  PSD  for  wavenumbers  greater  than  the 
inverse  Fresnel  distance  to  a much  flatter  behavior  for  larger  o^. 
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Wavenumber  (km~^) 

?1 . Normalized  intensity  PSD  for  exponential  phase-screei 
= 3 km.  A = 0.04  m. 


4.3 


Simulation  Results  - 10  kilometer  Correlation  Length 


I'or  phase  correlation  length  = 10  km,  the  two  phase-screen 
PSD's  generated  by  the  code  differ  significantly  in  the  wavenumber  region 
near  the  inverse  Fresnel  length  as  shown  in  Figure  22.  This  difference  is 
dramatically  shown  in  Figures  23  and  24  which  show  realizations  of  intensity 
for  0^  = 10.0  and  31.6  radians  for  the  power-law  PSD  and  for  = 10.0, 
50.0,  and  100  radians  for  the  exponential  PSD.  The  corresponding  phase 
realizations  are  shown  in  Figures  25  and  26. 

Figures  27  and  28  show  the  moan  intensity  PSD  corresponding  to 
intensities  shown  in  Figures  23  and  24. 
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Figure  22.  Comparison  of  simulation  generated  phase-screen  PSO's  for 
correlation  length  ^ = 0.04  m.  (o^.  = 10 

radians  in  figure).  ^ 
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Figure  26.  Realizations  of  signal  phase  for  exponential 
PSD  with  2.^  = 10  km. 
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Figure  27.  Normalized  intensity  PSD  for  power-law  phase-screen  PSD 
with  «,  = 10  km,  A = 0.04  m. 
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SECTION  5 


5.0  CONCLUSIONS 

This  report  has  presented  and  analyzed  simulated  results  for 
electric  field  propagation  through  a thick,  extended  region  of  high-altitude 
striations  of  electron  density.  The  simulated  results  correspond  to  the 
pi'oblem  of  communication  from  a DSCS  satellite  to  a ground  station  at  a 
frequency  of  7.5  GHz.  The  results  arc  parameterized  in  terms  of  the  phase- 
screen  power  spectral  density  (PSD)  and  phase  standard  deviation  which 
may  be  related  to  the  in-situ  electron  density  PSD  and  electron  density 
fluctuation  0^^  by  formulae  developed  in  Appendix  B.  Statistical  results 
are  presented  for  S^j  scintillation  index  and  intensity  correlation  length 
for  a power- law  and  an  exponential  phase-screen  PSD  for  values  of  phase 
standard  deviation  ranging  from  1 to  100  radians  and  phase  correlation 
length  ranging  from  1 to  10  km. 

Examples  of  specific  realizations  of  the  electric  field  amplitude 
and  phase  at  the  receiver  are  also  sho\\rn.  These  realizations  of  signal 
structure  are  being  used  in  simulations  of  DSCS  receiver  operation,  reported 
elsewhere. 
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APPENDIX  A 


MULTIPLE  PHASE  SCREEN  CODE 


This  Appendix  describes  the  multiple  phase  screen  propajjation  code 
as  used  to  solve  the  scalar  wave  equation.  Also  included  is  a description 
of  the  technique  used  to  generate  the  random  phase-screens. 

A.l  Scalar  Wave  Equation  Solution 

The  problem  of  interest  is  that  of  a plane  wave,  incident  at 
z = 0 on  a random  medium  consisting  of  index  of  refraction  fluctuations. 
Assiune  that  the  random  medium  extends  from  z = 0 to  some  distance  and  that 
its  autocorrelation  function  is  not  a function  of  y.  That  is,  the  medium 
can  be  thought  to  consist  of  a number  of  random  striations  or  tubes  which 
are  infinitely  extended  in  the  y-direction. 


The  propagating  wave  is  then  only  a function  of  x and  z. 
Ignoring  depolarization,  the  scalar  wave  equation  nmy  bo  written 


2 2 2 

(V  + k n )i|»  = 0 


(A-l) 


where 


and  n = 1 + An  is  the  index  of  refraction  with  An  a small  perturbation. 
The  equation 


^ ^ ^ 2k^Anu  = 0 (A-2) 

3Z  ^ 3x 
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may  be  obtained  from  Equation  A-1  using  the  substitution 


i|j  = u expCikz) 


(A-3) 


and  ignoring  terms  in  (An)“. 

The  complex  amplitude  u can  vary  markedly  over  distances  no 
smaller  than  the  inhoiiiogenoity  scale  size  Jt;  that  is,  the  worst  case 

variation  of  u in  tlie  direction  of  propagation  can  be  characterized  by 

2 2 2 2 
u ~ cxp(-z  wlioro  2:  Jl.  Thus,  the  second  derivative  3 u/3z  is 

of  the  order  of  On  the  other  hand,  the  second  term  of  Equation  A-2 

is  of  the  order  of  u/Xl  . Thorefoi'c,  for  A <<  £ , or  A <<  £,  the  first 

m 111 

term  of  Equation  A-2  may  be  ignored  with  I'espect  to  the  second  term  to 
obtain  the  parabolic  wave  equation 


2 

^ + i2k  ~i  + 2k^Anu  = 0 
3x“ 


(A-4) 


Now  for  a given  initial  field,  the  px'opagation  simulation  must 
generate  values  of  the  field  such  that  Equation  A-4  is  satisfied.  To 
do  this,  the  medium  is  divided  into  screens  defined  by  planes  on  which  z 
is  constant.  To  propagate  from  one  scx’een  to  the  next  write  u(x,z)  in 
the  form 


u = we 


(A-5) 


wJiere 


so  that 


z+Az/2 

<p(x,z)  = ik  J'  dz'  An(x,z') 
z-Az/2 


(A-6) 


9(})  _ 


ikAn(x,z) 


CA-7) 
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Equations  A-5  and  A-7  may  be  substituted  into  Equation  A-4  to  obtain 


(we*^)  + 2ik  = 0 


(A-8) 


To  solve  this  equation,  assume  that  the  step  size  Az  is  sufficiently 
small  so  that  the  effect  of  the  exponential  factors  in  Equation  A-8  is 
small.  Thus,  one  can  solve 


9 w 9w  . 

^ ♦ 2ik  = 0 

)X 

A simple  solution  may  bo  obtained  using  the  Eoui'icr  transfoi'in  relation 

-ikx 


(A-9) 


k(x,:) 


J dK  W{K,z 


)o 


(A-10) 


which  may  bo  substituted  into  Equation  A-9  to  obtain 


A 


(A-11) 


whose  solution  is 


W(K,z) 


1S'{K,0) 


or  using  Equation  A- 12  in  Equation  A-10 


w(x,zl  = y'dK  W(K,0)  exp  |-  - i^xj 


(A-12) 


(A-13) 


Equation  A-13  may  be  rewritten  by  applying  Equation  A-12  for  z = z^  and 
z = where  Zj  and  are  two  successive  phase -screen  locations 
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00  2 
f ^ / iK  (Zo-zJ 

= /dKK(K,zp  exp  I 2k 


(A-14) 


Now,  using  Equation  A-5  and  Equation  A-3  and  replacing  by  the  electric 
field  we  obtain  from  Equation  A-14 

ikCZp-z.)  r A.  / j \ 

E(x,Z2)  = e / dK  E(K,zp(exp  - (z2-zp  - iKxJ  (A-15) 


whore  E(K,zp  is  now  the  Fourier  transform  of  the  product 


z,+Az/2 


E(x,zp  exp 


dz'  An(x,z')] 


(A-16) 


Zj-Az/2 


At  this  point  it  is  interesting  to  note  that  one  can  obtain  the 
Fresnel-Kirchhoff  integral  by  substituting  the  Fourier  transform 


E(K 


d?  E(C,z^)o 


iKC 


into  Equation  A-15  and  perfoi'ming  the  K integration  (Gradshtoyn  and 
Ryzhik,  p.  397)  as 

•1/2  00 


h(x,Z2) 


^i2Ti(z2-Zj)  J j 


dCe 


ik(z2-Zj) 


exp  j i -j  k(x-C)^/(z2-Zj)  j n(^,zp  (A-17) 


The  propagation  simulation  thus  breaks  up  the  random  medium  into 
a number  of  screens.  The  electric  field  incident  on  a phase-screen  in  then 
multiplied  by  the  phase  shift  caused  by  the  random  medium  lying  between  the 
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previous  and  current  phase-screen  (Equation  A-12).  The  wave  then  propagates 

as  if  in  free  space  to  the  next  phase-screen  (Equation  A-11).  Propagation 

from  one  phase-screen  to  the  next  is  accomplished  very  efficiently  by  the 

use  of  the  fast  Fourier  transform  so  that  once  the  electric  field  is 

multiplied  by  the  appropriate  complex  phase  term  (as  in  Equation  A-12)  its 

2 

Fourier  transform  is  taken,  multiplied  by  exp(-iK  and  then 

another  fast  Fourier  transform  is  performed  to  obtain  the  electric  field 
at  the  next  screen  E(x,Z2). 

A. 2 Generation  of  the  Random  Phase  Screen 

The  following  procedure  taken  from  Buckley  (1975)  was  used  to 
generate  the  random  phase  screen  with  arbitrary  variance  and  power  spectral 
density. 


Assume  that  P^  is  a sequence  of  random  numbers  uniformly  distri- 
buted between  -1/2  and  1/2.  Then 

(P.)=  0 

("ii’j)  ■ h hi 

where  is  the  Kronecker  delta  function.  Now  choose  the  phase  as  a 

sequence  given  by  the  circular  (see  Enslein,  et  al.,  p.  380)  correlation 

N-1 

4>.  = <p(ZAx)  = T w P ; Z = 0,  N-1  (A-19) 

^ j=0  ^ J 

where  the  phase  screen  is  sampled  at  increments  Ax  = L/N  where  L is  the 
total  length  of  the  screen.  The  w^  are  a set  of  weights  which  are  to  be 
determined. 
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Since  each  value  (j)j^  is  a sum  of  a large  number  N of  independent 
random  variables,  it  follows  from  the  central  limit  theorem  that  the  sequence 
{(p^)  will  possess  a Gaussian  distribution  with  zero  mean.  Now  the  phase 
autocorrelation  function  Pj,  may  be  formed  by  noting 


N-1  N-1 

= "'3  ^3+i 


'^2,  ‘^2,+i+k 


) 


N-1 

= I, 

1 


0 


1 oT 


N-1 

E 


£=0 


"'z  12  '^j,2-+k 


^ N-1 

= 12  g,  Vk 


(A-20) 


Set 


to  obtain  the  autocorrelation  function 


(A-21) 


_ 1 

■ .2 


1 

Pv  = 9 E '’'2.  "'2,+k 

12({.^  Z%  ^ 

The  p.,  sequence  may  be  interpreted  as 


k=0,...N-l 


(A-22) 


p(kAx),  k = 0,...N-1;  that  is,  Pj, 


is  a discrete  representation  of  the  continuous  autocorrelation  function. 


I 


i 
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Now  in  order  to  relate  the  still  unknown 
to  the  power  spectral  density  of  the  phase  we  must 
pair  00 

p(x)  = y's(K)e"^‘^^dK 
-00 


weights  {Wj^,  k=0,...N-l} 
use  the  Fourier  transform 

(A-23) 


SCK) 


(A-.24) 


and  their  discrete  representations 


p(kAx) 


N-1 

2TrAQ  Y,  sa2TTAQ) 
£=0 


-i(kAx)(£2TiAQ) 

® * 


k=0,...,N-l 


(A-25) 


S()l2TrAQ) 


k=0 


A=0,..,N-1 


(A- 26) 


whore  AK  = 2TrAQ.  Noting  that  Ax  = L/N,  where  L is  the  length  of  a 
phase  screen,  AQ  = 1/L,  AxAQ  = 1/N,  Equation  A-25  and  Equation  A-26  may  be 
rewritten  as 

p.  = p(kAx)  = 27rAQ  S,  k=0,...N-l  (A-27) 

^ il=0  ^ 

= 6(£27rAQ)  = Pj.  £=0,...N-1  (A- 28) 

k=0 


Now  the  Fourier  transfer  of  Equation  A-22  is  given  by 


S 

r 


Ax 

2n 


124. 


N-1  N-1 

E E 

k=0  1=0 


w > 


il+k 


i27rkr/N 

e 


(A-29) 
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The  Fourier  transform  of  the  weight  is  given  by  IV  as 


Wfl  = 2ttAQ  V W e 
X.  ^ m 

m=0 


-i2vmi/N 


(A-30) 


and  similarly  for 


N-1 


«Uk  = 2''4Q  »„  = 

n=0 


■i2Tr(Uk)n/N' 


(A-31) 


Now  since  the  weights  {Wj,}  are  real  complex 

conjugate  of  Equation  A-30  and  Equation  A-31  with  Equation  A-29  to  obtain 


2,, ^,2 


r 2 ^ ^ ^ m 

2‘\v^q  k i m n 


^-i2Tr()l+k)n/N  ^i27ikr/N 
' n 


()*  +i27rmii/N 

Iv  ft 


(A-32) 


Now  since 


-27r)ln/N  i2TrmVN 


2^ 

£=0 


N6 


inn 


N-1 


-i2Trkn/N  i2'nkr/N 

ft  e 


T c C = No 


k=0 


nr 


Equation  A-32  may  be  summed  to  yield 


s = nAQ)-  Ax  .r  j-  I 


(A-33) 


(A- 34) 


or  since  NAQ  = 1/Ax  the  Fourier  transform  of  the  weights  is  given  by: 


^ /bAx({)  “ 

IV  = / L.  s 

^ Tf  r 


(A-35J 


where  S is  the  discrete  Fourier  transform  representation  of  the  phase- 
s-'.reen  power  spectral  density. 
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Thus,  to  generate  the  phase  values  for  a phase-screen,  one  may 
start  with  the  power  spectral  density,  obtain  the  Fourier  transform  of  the 
weights  via  Equation  A-35,  and  take  the  inverse  Fourier  transform  to  obtain 
the  weights  themselves. 
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APPENDIX  B 


RELATIONSHIP  BETWEEN  PHASE  AND  ELECTRON  NUMBER  DENSITY 


B.l  Phase  Screen  PSD  Related  to  Electron  Density  PSD 

As  discussed  in  Appendix  A,  the  electric  field  at  each  phase 
screen  is  multiplied  by  the  factor  where 


(j)(x)  = L / An(x,z)  dz 

-Lz 

2 


2 


(B-1) 


which  effectively  collapses  the  index-of-refraction  fluctuations  in  the 
region  from  -Az/2  to  Az/2,  In  the  absence  of  significant  electron 
collisions,  the  index  of  refraction  is  given  as 


2 1 
n = 1 


N 


_e 

★ 


(B-2) 


where  is  the  electron  number  density  and  n is  the  critical  electron 
number  density 


* C-m  OJ  R 9 X 

n = = 1.24x10"“  f [eVcm'^l 


(B-3) 


where  is  the  permittivity  of  free  space,  m^  is  the  mass  of  an  electron, 

e is  the  electron  charge,  and  f is  the  radio  frequency  in  Hertz.  From 

* 

Equation  B-2,  the  deviation  in  refractive  index  is,  for  N « n 

© 


, AN 

A 1 e 

An  = 2 * 

n 


(B-4) 
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PRSCEDIKG  PAGE  BU’.<K-NOT  FIUffiD 


so  that  Equation  B-1  can  be  written  as 

Az 
2 


j f AN^(x,z)  dz 

" -Az 


(B-5) 


Now  the  autocorrelation  function  can  be  formed  as 


aj  = (A(()(x)  A4>(x+5)) 


Az  Az 

T T 


f f (ANg(x,z}  An^(x+5,z')  ) dzdz' 
n •<. 


-Az  -Az 
2 2 


~ „ Az  Az 

Om  2 2 

-T^// 

n A *'a 

-Az  -Az 
2 2 


B^  (5,z-z')  dzdz' 


(B-6) 


The  double  integral  here  can  be  reduced  to  a single  integration  by  a change 
of  variables  as  explained  by  Pan  ■ ’.is  (1965),  page  325  as 


2,2  . 

0.,  k Az 


4 n ' e 

-Az 


(B-7) 


If  Az  is  greater  than  the  co.'relation  xength  of  the  electron  density 
fluctuations,  then  the  contribution  of  the  second  term  in  the  integration 
will  always  be  negligible  so  that 
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i = 


Az 


oo 


(5,n)  dn 


(B-8) 


iK  5 

Now,  by  multiplying  both  sides  of  Equation  B-8  by  e /In  and  integrating 
with  respect  to  C and  using  the  Fourier  transform  relationship  given  by 
Equations  B-11  and  B-14  and  their  2-D  extensions  one  obtains 

♦ (K  ) - = « 

n e 


Equation  B-9  is  a general  relationship  between  the  phase  spectrum  and  the 
2-dimensional  electron  density  spectrum. 


Setting  f = 0 in  Equation  B-8  and  noting  that  B^(0)  = 1 gives 
the  relationship  between  the  variance  of  phase  and  the  variance  of  electron 
density  fluctuations 


2 


a 


4' 


2 

^N 


(o,n)  dn 


(B-IO) 


B.2  One,  Two  and  Three  Dimensional  Power  Spectra 


Given  an  isotropic  and  homogeneous  random  medium  with  the  auto- 
correlation functions  Bj(x),  B^(p),  and  B^Cr)  in  one,  two,  and  three- 
dimensions,  then  the  respective  power  spectral  densities  are: 


^ (K)  = Ifj-  J dx  cos  Kx  B^(x) 


4>2(K) 


2 r 

^ J P dp  Jq(Kp)  B2(p) 


(B-li) 


(B-12) 
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2 

(f  (K)  = / dr  r sin  Kr  B (r) 

27r  K ^ 


(B-13) 


which  have  the  inverse  transforms: 


/« 


0 Bj(x)  = J dK  cos  Kx  «>j(K) 


(B-14) 


OO 

'/ 


O^-B^Cp)  = 271  j dK  K Jq(Kp)  4>2(K) 


(B-15) 


o^BjCr)  = " f clK  K sin  Kr  <I>j(K) 


(D-16) 


where  the  subscripts  again  deiiote  the  dimensionality. 

Those  relationships  above  are  consistent  with  the  relationships 
between  the  specti'a  in  various  dimensions  (Rufenach,  1975): 


4>^(K  ,K  ) = I dK  <I>,(K  ,K  ,K  ) 
2^  x’  y*  J z 3^  X’  y’ 


(B-17) 


4>,  (K  ) = / dK  4)-,(K  ,K  ) 

y x-*  J y 2'-  x’  y* 


(B-18) 
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Under  the  assumption  that  the  autocorrelation  function  has  the 
same  form  in  one,  two,  or  three  dimensions,  the  following  relationships 
are  easily  shown  for  some  common  autocorrelation  functions.  Only  the 
3-dimensional  form  of  the  autocorrelation  function  is  specified. 


1. 


Exponential  Autocorrelation  Function 


B(r)  = exp  j -Irl/iQ j 
“n  r 

‘ ir-TTr 


1+K  r: 


0 


of,  2 

N 1' 

(j)  /•[;■)  s 2.  Q 

' (1^k2  rj)  V2 

cl 

, ...  ^ ^0 

“ ~2'  ,..,,2  2,2 


^ (1+r  rj)‘ 


(B-IS) 


(B-20) 


(B-21) 


(U-22) 


2.  Power-Law  Autocorrelation  Functions 


B(r)  = 


(B-23) 
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(n  need  not  be  an  integer) 


4-n 

2^  r 


= d 


n-4 

2 


® 7r^/2  J,(-^-)  2 


(B-24) 


3-n 


^3(K)  = oj 


2 2 

r' 

c 

TT  r ( 

2 / 

2-n 

0 2 

■3 

2 

r 

0 

IW 

TT 

r(Jj 

n-5 

2 


nr-S 

~2 


(Kro), 


n-6 

('<'■0)  ' ''n.6 


(B-25) 


(B-26) 


where  K(x)  is  the  Bessel  function  of  the  third  kind, 

_9 

B.3  Two  Dimensional  K “ in-situ  Spectrum 

(K)>  the  one-dimensional  power  law  phase  power  spectral  density 
’>sud  in  this  report  may  be  derived  by  assuming  an  in-situ  K ^ spectrum  for 
electron  density  fluctuations  as  given  by  Equation  B-20.  Then,  Equation  B-21 
and  Equation  B-9  yield  the  phase  power  spectral  density 


N. 


0^(K)  = 


2TrAz  k ^ 

^ *2  2n  „2  2,  3/2 
4n  (1+K  r ) 


(B-27) 


The  integral  specified  by  Equation  B-10  may  be  performed  using  Equation  B-19 
as  the  autocorrelation  function  to  obtain 


2 '’0  2 

% " - *2 
^ 2n  e 


(B-28) 
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or  using  Equation  B-28  in  Equation  B-27  one  comes  full  circle  to  obtain 


f ^7:^  3/2 


which  is  identical  to  Equation  1-1  with  r^  replaced  by  Lq, 


B.4  Two-dimensional  Exponential  in-si tu  Spectrum 


121 

(K),  the  one-dimensional  exponential  phase  PSD  given  by  Equa- 
tion 1-2  may  be  obtained  from  Equation  B-25  with  n = 6 so  that  Equation 
B-25  represents  a two-dimensional  exponential  PSD  for  electron  density 
fluctuations 


e 


2 

2 

2ii  ® 
e 


The  autocorrelation  function  is  then 


B,  (rl  = 


e (l  + Cr/rJ") 


2,  5/2 


(B-30) 


(B-31 


So  that  the  phase  PSD  may  be  determined  as 


where 


2 It  Az  2 

' , *2  ‘’n 

^ 2n  e 
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(B-32) 


(B-33) 


APPENDIX  C 


PHASE-SCREEN  REALIZATIONS 


For  all  the  7.5  Gilz  DSCS  simulation  results  presented  in  this 
report  the  striated  region  was  modeled  by  ton  phase-screens  and  the  simu- 
lation was  run  ten  times,  each  time  with  a different  set  of  phase-screen 
realizations,  to  obtain  statistical  results,  rb.is  Appendix  shows  all  ten 
phase-screens  used  for  the  first  realization  and  also  the  observed  intensity 
at  each  screen  so  that  the  reader  may  "watch"  as  the  wave  propagates  through 
the  medium  to  the  receiver  plane.  For  this  case  the  exponential  PSD  given 
by  Equation  1-2  with  a phase  standai'd  deviation  = 31.6  radians  and  a 
correlation  length  = 3.0  km  is  used. 

The  10  phase  screens  arc  shown  in  Figure  C-1  (a-d)  as  a function 
of  distance  z from  the  top  of  the  disturbed  region  (see  Figure  1 (a)). 


Figure  C-2  (a-d)  shows  the  intensity  observed  at  the  second 
through  tenth  phase  screen  locati  i (the  intensity  at  the  first  phase- 
screen  location  is  a constant)  anu  at  the  receiver  plane  located  300  km 
from  the  tenth  phase  screen.  It  is  interesting  to  try  to  trace  the  focus- 
ing effects  which  are  occurring  as  the  wave  propagates  through  the  screen. 
The  Gaussian  irregularity  denoted  by  an  arrow  on  the  first  phe-se  screen 
(z  = 0,  Figure  C-1  (a))  can  be  seen  to  be  the  cause  of  the  intensity  dif- 
fraction pattern  which  is  traceable  for  z = 3,300  km  to  z = 15000  km 
and  denoted  by  arrows  on  each  observation  screen. 


This  realization  of  all  ten  screens  can  be  used  as  another  visual 


I 


} ' 
? I 


test  of  the  accuracy  of  the  simulation  by  using  the  Fresnel-Kirchhoff 
equation  to  analytically  compute  the  diffraction  pattern  of  the  Gaussian 
lens  shown  in  Figure  C-3.  The  diffraction  pattern  is  shown  as  a function 
of  z on  the  same  scale  as  used  in  Figure  C-2  and  shows  a striking 
resemblance  to  the  results  generated  for  the  random  screen  realization. 


Phase  (radians) 


Figure  C-l(d).  Realizations  of  phase-screens  with  exponential  PSD 
(o,  = 31 .6  rad,  = 3 km). 


{ 


Figure  C-2  (d).  Realizations  of  intensity  at  each  observation  screen 
location. 
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Figure  C-3.  Analytic  diffraction  pattern  of  a Gaussian  lens  as  a function 
of  distance  Z from  the  lens. 
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Director 

TRASANA 

ATTN:  ATAA-SA 

ATTN:  TCC,  F.  Payan,  Jr. 

ATTN:  ATAA-TAC,  LTC  John  Hesse 

Director 

U.S.  Array  Ballistic  Research  Labs. 

ATTN:  Lawrence  J.  Puckett 

Commander 

U.S.  Array  Corrm-Elcc.  Engrg.  Instal.  Agy. 
ATTH:  EED-PEO,  George  Lane 

Commender 

U.S.  Arra/  Electronics  Coirmand 

ATfH:  DRSEL-NI.-RD,  H,  S.  Bcnnet 
ATTH:  CR$EL-PL-EHV,  Hans  A.  Bomke 

Cortmander 

U.S.  Array  Foreign  Science  & Tech.  Ctr. 

ATTN:  R.  Jones 
ATTH:  P.  A.  Crowley 

Commander 

U.S.  Array  Missile  Command 

ATTN:  DRSMI-YTT,  W.  G.  Preussel,  Jr. 
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DEPARTMENT  OF  THE  AIR  fORCE 


Commander 

U.S.  Artny  Missile  Intelligence  Agency 
ATTN:  Jim  Gamble 

Conmander 

U.S.  Am\y  Nuclear  Agency 

ATTN:  MONA-WE,  0.  Berberet 

Commander 

U.S.  Anny  Satcom  Agency 
ATTN:  Doc.  Con. 

DEPARTMENT  OF  THE  NAVY 

Chief  of  Naval  Operations 

ATTN:  OP  943,  LCDR  Huff 
ATTN:  Ronald  E.  Pitkin 
ATTN:  Alexander  Brandt 

Chief  of  Naval  Research 
ATTN:  Code  418 
ATTN:  Code  461 

Commander 

Naval  Electronic  Systems  Cotimand 

Naval  Electronic  Systems  Cmd.  Hqs. 

ATTN:  PME  IIT-T,  Satellite  Comm.  ProJ.  Off. 
ATTN:  NAVAL EX  034.  T.  Barry  Hughes 
ATTN:  PME  117 

Connanding  Officer 

Naval  Intelligence  Support  Ctr. 

ATTN:  STIC  12.  Mr.  Dubbin 

Conmander 

Naval  Ocean  Systems  Center 

ATThi:  Code  0230,  C.  Baggett 
ATTN:  Code  532.  Hilliam  F.  Holer 

3 cy  A'  .;  Code  2200 
AlIN:  R.  Eastman 

Director 

Naval  Research  Laboratory 

ATTN:  Code  7701,  Jack  0.  Brown 
ATTN:  Code  5400.  Hg.  Conm.  Dir..  Brace  Hald 
ATTN:  Code  7700,  Timothy  P.  Coffey 
ATTN:  Code  5460,  Electromag.  Prop.  Br. 

ATTN:  Code  5461,  Trans.  Inon.  Prop. 

ATTN:  Code  5430,  Satellite  Comm. 

ATTN:  Code  9465,  Prop.  Applications 

Commander 

Naval  Space  SurveT  lance  System 
ATTN:  CAPT  J.  h.  Burton 

Officer- in-Charge 

Naval  Surface  Weapons  Center 

ATTN:  Code  WA501,  Navy  Nuc.  '’rgrns.  Off. 

Conmander 

Naval  Surface  Weapons  Center 

ATTN:  R.  F.  Butler,  OF- 14 

Director 

Strategic  Systems  Project  Office 

ATTN:  NSSP-2722,  Fred  Wimberly 
ATTN:  •iSP-2141 


Commander 

ADC/DC 

ATTN:  DC,  Mr.  Long 

Commander 

ADCOM/XPD 

ATTN:  XPQDq 

AF  Geophysics  Laboratory,  AFSC 
ATTN:  PHD,  John  P.  Mullen 

ATTN:  OPR,  Alva  T.  Stair 

ATTN:  OPR,  James  C.  Ulwick 
ATTN:  LKB,  Kenneth  S.  W.  Champion 
ATTN:  PHD,  Jurgen  Buchau 
ATTN:  OPR,  Harold  Gardner 
ATTN-  SUOL,  Rsch.  Lib. 

AIi.N:  PHP,  Jules  Aarons 

AF  Weapons  Laboratory,  AFSC 
ATTN:  SUL 

ATTN:  CA,  Arthur  H.  Guenther 
ATTN:  SAS,  John  M.  Kamra 

ATTN:  DYC.  Capt  L.  Hittwer 

AFTAC 

ATTN:  TF,  Maj  Wiley 
ATTN:  TN 

Air  Force  Avionics  Laboratory,  AFSC 
ATTN:  AAD,  Allen  Johnson 

Headquarters 

Electronic  Systems  Divislon/XR 
ATTN:  XRC,  Lt  Col  J.  Morin 
ATTN:  XRE,  Lt  Michaels 

Headquarters 

Electronic  Systems  Division/YS 
ATTN:  YSEV 

Headquarters 

Electronic  Systems  Division,  AFSC 
ATTN:  Jim  Oeas 

Commander 

Foreign  Technology  Division,  AFSC 
Anti:  NICO  Library 
ATTN:  ETD,  B.  L.  Ballard 

Hq.  USAF/RD 

ATTN:  RDQ 

Commander 

Rome  Air  Development  Center,  AFSC 
ATTN:  V.  Coyne  Ocse 
ATTN:  EHTLO.  Doc.  Lib. 

Commander 

Rome  Air  Development  Center,  AFSC 
ATTN:  ETEI,  A.  Lorentzen 

SAHSO/HN 

ATTN:  HNNL,  Lt  Col  Kennedy 
SAHSO/SK 

ATTN:  SKA,  Lt.  Maria  A.  Clavin 
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DEPARTIIEHT  OF  THE  AIR  FORCE  (Continued) 
SAHSO/SZ 

ATTN:  S2J,  Major  Lawrence  Doan 
SAHSO/VA 

ATTH:  YAT,  Capt  L.  Blackwelder 

Commander  in  Chief 
Atrategic  Air  Command 
ATTN:  NRT 

ATTN:  XPFS,  Haj  Brian  G.  Stephan 
DEPARTMENT  OF  ENERGY 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Analytical  Systems  Engineering  Corp. 

ATTN:  Radio  Sciences 

The  Boeing  Company 

ATTN:  0.  Murray 
ATTH:  Glen  Keister 

Brown  Engineering  Company,  Inc. 

ATTN:  Romeo  A.  Deliberis 

University  of  California  at  San  Diego 
Marine  Physical  Lab.  of  the  Scripps 
Institute  of  Oceanography 


EG4G.  Inc. 

ATTN: 

Henry  G.  Booker 

Los  Alamos  Division 

ATTN 

James  R.  Breedlove 

Charles  Stark  Draper  Laboratory,  Inc. 

ATTH: 

J.  P.  Gilmore,  MS  63 

University  of  California 

Lawrence  Livermore  Laboratory 

ATTH: 

D.  B.  Cox 

ATTN 

Ronald  L.  Ott,  L-531 

Computer  St • 

. ■ s Corporation 

ATTN 

Tech.  Info.  Dept. , L-3 

A"1N 

1 lank 

ATTN 

Donald  R.  Dunn,  L-156 

ATTN 

Ralph  S.  Hager,  L-31 

ESL,  Inc. 

i 

ATTN: 

J.  Roberts  , 

Los  Alamos  Scientific  Laboratory 

ATTN: 

C.  H.  Prettie  1 

ATTN 

Doc.  Con.  for  John  Zinn 

ATTN: 

James  Marshall 

ATTN 

Doc.  Con.  for  R.  F.  Taschek 

ATTN: 

V.  L.  itower  i 

ATTN 

Doc.  Con.  for  John  S.  Malik 

ATTN 

Doc.  Con.  for  Eric  Jones 

Ford  Aerospace  & Conriunl cat ions  Corp.  ( 

ATTN: 

J.  T.  Hattingley,  MS  X22  j 

Sandia  Laboratories 


ATTN: 

Doc. 

Con. 

for 

0. 

ATTN: 

Doc. 

Con. 

for 

J. 

ATTN: 

Doc. 

Con. 

for 

T. 

ATTH: 

Doc. 

Con. 

for 

A. 

Org. 

1245 

ATTH: 

Doc. 

Con. 

for 

W. 

OTHER  GOVERNMENT  AGENCIES 

Department  of  Cormerce 
Office  of  Telecommunications 
Institute  for  Telecom.  Science 


ATTH 

ATTN 

ATTN 


L.  A.  Berry 
G.  Reed 

William  F.  Utlaut 


NASA 

Godoard  Space  Flight  Center 

ATTN:  ATS-6,  OFC  P.  Corrigan 

National  Oceanic  i Atmospheric  Admin. 
Environmental  Research  Laboratories 
Oepartm.ent  of  Fommerce 


ATTH 

ATTN 

ATTN 


Joseph  H.  Pope 
C.  L.  Rufenach 
Richard  Grubb 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 
Aerospace  Corporation 


AT  IN 
ATTH 
ATTH 
ATTN 
ATTN 
ATTN 
ATIN 
ATTN 
ATTN 


'j.  P.  Bower 
T.  H.  Salmi 
Norman  0.  Stockwell 
F.  A.  Horse,  A6,  Rm.  2407 
V.  Josephson 

J.  E.  Carter,  120,  Rm.  2209 
Irving  H.  Garfunke'. 

F.  E.  Bond,  Al,  Rm.  5003 
0.  P.  Olsen.  120,  Rm.  2224E 


General  Electric  Company 

Space  Division 

Valley  Forge  Space  Center 

ATTN:  M.  H.  Bortner,  Space  Sci.  Lab. 

General  Electric  Company 

TEMPO-Center  fo.'  Advanced  Studies 
ATT‘i.  Tir.i  Stephens 
ATTN:  Tom  Barrett 
ATTH:  B.  Gambill 
ATTH:  William  McNamera 
ATTN:  Warren  S.  Knapp 
ATTN:  Hack  Stanton 
ATTN:  Don  Chandler 
ATTN:  DASIAC 

General  Electric  Company 
ATTN:  F.  A.  Reibert 

General  Research  Corporation 
ATTH:  Joel  Garbarino 
ATTH:  John  Ise,  Jr. 

Geophysical  Institute 

University  of  Alaska 

ATTN:  T.  H.  Davis 
ATTH:  Heal  Brown 
ATTN:  Tech.  Lib. 

GTE  Sylvania,  Inc. 

Electronics  Systems  Grp.,  Eastern  Div. 
ATTH:  Marshal  Cross 

HSS,  Inc. 

ATTH:  Donald  Hansen 

Inti.  Tel.  4 Telegraph  Corporation 
ATTN:  Tech.  Lib. 
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DEPARTKliHT  OF  DEFENSE  CONTRACTORS  (Continued) 

Institute  for  Defense  Analyses 
ATTN:  Ernest  Bauer 
ATTN:  Joel  Bengston 
ATTN:  Hans  Wolfhard 
ATTN:  J.  (1.  Aein 


ATTN:  S.  R.  Goldman 

Johns  Hopkins  University 

Applied  Physics  Laboratory 

ATTN:  Document  Librarian 
ATTN:  ThOMS  Potcmra 
ATTN:  John  Dassoulas 

Kaman  Sciences  Corporation 
ATTN:  B.  J.  Bittner 
A^tn:  T.  Meagher 

Linkabit  Corporation 

ATTN:  Irwin  Jacobs 

Lockheed  Missiles  A Space  Co.,  Inc. 

ATTN:  D.  R.  Churchill 
ATTN:  Dept.  60-12 

Lockheed  Missiles  4 Space  Co.,  Inc. 

ATTN:  Richard  G.  Johnson,  Dept.  52-12 
ATTN:  Martin  Halt,  Dept.  52-10 

M.I.T.  Lincoln  Laboratory 

ATTN:  Mr.  Walden,  XU3 
ATTN:  James  H.  Pannell,  L-246 
ATTN:  David  M.  Towle,  Lib, , A-082 
ATTN:  0.  Clark 

McDonnell  Douglas  Corporation 
ATTN:  J.  Houle 
ATTN:  George  Mrot 
ATTN:  William  Olson 
ATTN:  N.  Harris 

Mission  Research  Corporation 

ATTN:  Warren  A.  Schlueter 
ATTN:  Raloh  Kilb 
ATTN:  R.  Hendrick 
ATTN:  R.  Bogusch 
ATTN:  Dave  Sowle 
ATTN:  F.  Fajen 
ATTN:  H.  F.  Crevicr 
ATTN:  Steven  L.  Gutsche 
ATTN:  D.  Sappenfield 
ATTN:  0.  L.  Knepp 
5 cy  ATTN:  Tech.  Lib. 


The  Mitre  Corporation 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Pacific-Sierra  Research  Corp. 

ATTN:  E.  C.  Field,  Jr. 

.otometrics,  Inc. 

ATTN:  Irving  L.  Kofsky 

Physical  Dynamics,  Inc. 

ATTN:  Joseph  B.  Workman 
ATTN:  A.  Thompson 

R 4 D Associates 

ATTN:  William  J.  Karzas 
ATTN:  William  B.  Wright,  Jr. 

ATTN:  Forrest  Gilmore 
ATTN:  Robert  E.  Lelevier 
ATTN:  Bryan  Gabbard 

The  Rand  Corporation 

ATTN:  Ed  Bedrozian 
ATTN:  Cullen  Crain 

Raytheon  Company 

ATTN:  Barbara  Adams 

Science  Applications,  Inc. 

ATTN:  E.  A.  Straker 
ATTN:  Daniel  A.  Hamlin 
ATTN:  Curtis  A.  Smith 
ATTN:  D.  Sachs 
ATTN:  Lewis  M.  LInson 
ATTN:  Jack  McOougall 

bcicncc  Applications,  Inc. 

Huntsville  Division 

ATTN:  Dale  H.  OiviS 

SRI  International 

ATTN:  Alan  Burns 
ATTN:  Halter  G.  Chesnut 
ATTN:  Charles  L.  Rino 
ATTN:  Ray  L.  Leadabrand 
ATTN:  Halter  Jayc 
ATTN:  H.  Baron 
ATTN:  L.  L.  Cobb 
ATTN:  Donald  Neilson 

Technology  International  Cor.noration 
ATTN:  W.  P.  Boquist 

TRW  Defense  4 Space  Sys.  Group 

ATTN:  R.  K.  Plebuch,  Rl-2078 

Visidyne,  Inc. 

ATTN:  Charles  Humphrey 
ATTN:  J.  H.  Carpenter 


C.  E.  Callahen 
Chief  Scienlivt , H.  Sen 
G.  Harding 
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